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I review recent lattice results on the spectrum and structure of baryons. Limita- 
tions due to the quenched approximation and un-physically heavy up and down 
quarks are discussed, and interfaces between first principles studies of QCD (or 
approximations to it) and "model building" are highlighted. 



Q ' 1. Introduction 

^3 ! 

pH \ With new exciting experiments ongoing and planned at several facilities, 

distributed over three continents, the past few years have witnessed the 
emergence of a new arena of research on the borderline between parti- 
' cle and nuclear physics: hadron physics. These developments have been 

■ complemented by several break-throughs in our theoretical understanding 

of the long distance, non-perturbative regime of hadronic matter, directly 
from QCD. Here I will focus on the most recent lattice results. 

Progress in the field has been slow but steady. Several "revolutions" in 
the lattice technology, for instance non-perturbative renormalisation and 
improvement programmes of action and operators, were often hardly iden- 
tifiable for the outside observer, a feature that we share with many other 
fields of research. In the case of strong QCD the failure of immediate sat- 
isfaction is due to the highly non-trivial phenomenology emerging from a 
Lagrangian as simple as that of QCD and yet so colourful and complex that 
even now it keeps experimentalists and theoreticians very busy. 

A major but all to often neglected application of lattice results is to 
establish validity range and applicability of phenomenological models or 
effective field theories directly from QCD and to calculate low energy pa- 
rameters that appear in these models or approximations. In quite a few sit- 
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uations, lattice predictions have matured to a degree that aUows direct con- 
frontation with experiment. Triggered by new experimental programmes 
and improved interaction with other theoreticians, an increasing number of 
lattice practitioners became interested in what has now become the domain 
of hadronic physics. This more focused approach resulted in an accelerated 
development of the field over the past two years with more to come in the 
future. I will go through the present state of lattice studies, discuss the 
two major sources of error, and review results on spectrum and structure 
of baryons, in particular A transition form factors, charge distribu- 

tions, three-body forces and the hyperon spin content. 

2. The state of the lattice 

Simulation results are obtained at finite lattice spacings and have to be 
extrapolated to the continuum limit. Moreover, the space-time box size 
should be big enough to comfortably accommodate the hadronic scales 
of interest. The dominant systematic uncertainties of present-day lattice 
results however tend to be related to unrealistically heavy u and d quark 
masses, which typically have to be chirally extrapolated over a long distance 
from > 0.25 mj, to the physical limit « 0.03 mj,, often combined with 
the quenched approximation, i.e. neglecting the polarisation effects of sea 
quarks on the QCD vacuum. 

The latter approximation, which easily reduces the computational effort 
by a factor of 10"^, is justified in the limit of the number of colours Nc ^ oo, 
however, it is not a priori clear whether 1/3 <C 1. Phenomenology^ suggests 
that this assumption is not as outrageous as it might sound at first. This 
seems substantiated by simulations of pure SU{Nc) gauge theories for^ 
A?'c = 2, . . . , 6 as well as by a comparison of quenched lattice predictions of 
the light hadron spectrum with experiment: ratios of light hadron masses 
on the lattice indeed have at last been found to be inconsistent with the 
observed spectrum."^ The differences are typically smaller than 10 %. 

However, there are cases in which "quenching" clearly matters. Quan- 
tum mechanical perturbation theory tells us that an O(e^) correction to 
an energy eigenvalue corresponds to an 0{e) correction on wave functions 
which hence can easily be changed by as much as 30 %. This will af- 
fect quantities like decay constants and parton distributions. The situation 
with for instance bottomonium fine structure splittings or clcctrowcak de- 
cay rates is even worse, ^ as these are roughly proportional to the square 
of the wave function at the origin. Strong decays do not exist in this ap- 
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Figure 1. Yearly SPIRES citations of Nielsen-Ninomya^ (NOGO) and GW** (GO). 



proximation as quarks and anti-quarks cannot annihilate, and therefore one 
might hesitate to expect very broad resonances to be modelled in a reason- 
able way. The same also holds true for any valence quark based model. We 
shall discuss some other shortcomings of quenching below. 

Fortunately, with several new supercomputer installations on the hori- 
zon by early 2004 we are at the brink of the fourth generation of simulations 
including sea quarks. These certainly will bring us down in the quark mass 
to about the point where the p meson becomes unstable and will at the 
same time allow for controlled continuum limit extrapolations. 

Recently a way has been found to implement the analogue of exact chiral 
symmetry without "doublers" on the lattice, known as overlap, Neuberger 
or Ginsparg- Wilson (GW) fermions.^ One particular (truncated) realisa- 
tion of the Neuberger operator are domain wall fermions.^ GW fermions 
are the only known consistent way of regularising a quantum field theory 
with massless fermions (like e.g. the standard model) and hence obviously a 
major theoretical break-through in particle physics as a whole. The famous 
Nielsen- Ninomiya no-go theorem^ is circumvented by relaxing the definition 
of chiral symmetry.^ I illustrate the yearly number of citations^ of these two 
important articles'''® published in 1981 and 1982, respectively, in Figure 1 



February 8, 2008 23:12 WSPC/Trim Size: 9in x 6in for Proceedings 



nstar 



4 

and leave judgement to the reader which paper would have better impacted 
on Research Assessment Exercises. In practical terms GW fermions mean 
that realistically light quark masses will become accessible to future numer- 
ical simulations, first qucnched"'^'' and subsequently un-qucnchcd. 

Going to light quark masses is expensive as the linear spatial lattice 
extent L has to be large relative to the inverse pion mass, to avoid strong 
finite size effects from pions exchanged "around" the boundaries. For in- 
stance if we intend to simulate a physical pion with say L > 4/m7r this 
implies L > 5.7 fm. The main limitation however is that the effort of in- 
verting the lattice Dirac operator scales with a large power of the inverse 
quark mass. This will improve with GW quarks that do not suffer from un- 
physical zero or near-zero modes: while at m,r > 0.5 GW fermions are 
easily two orders of magnitude more expensive to simulate than "standard" 
formulations, there will be a reward at small quark masses. In the quenched 
approximation this improvement factor has already turned out to be infi- 
nite, rendering a pion lighter than 200 MeV from impossible to possible. ^° 
Chiral fermions also greatly simplify operator mixing and renormalisation. 

3. Chiral extrapolations and "quenching" 

Based on the assumptions that QCD bound states are mesons and baryons, 
that there is a mass gap and spontaneous chiral symmetry breaking at zero 
quark mass, an effective low energy chiral Lagrangian can be derived in 
the spirit of the Born-Oppenheimer approximation. This Lagrangian will, 
to leading order, describe interactions between the (fast moving) massless 
Nambu-Goldstone pions and other hadrons. In nature quarks and thus 
pions are not massless and the leading mass corrections are formally of 
order m^/A-^sB where Ay^sB ~ 47rF^ > 1 GeV. Of course in an effective 
field theory the number of terms explodes at higher orders and predictive 
power is eventually lost, imless an early truncation is possible. 

Third generation lattice simulations with sea quarks were limited to un- 
realistically heavy pions, heavier than about 400 MeV. Only recently masses 
as low as 180 MeV have become possible in the quenched approximation,-'^'' 
albeit not yet at sufficiently large volumes. To allow for a first principles 
connection with phenomenology it is mandatory to establish an overlap 
region between lattice simulations and (leading order) chiral perturbation 
theory. In general the size of this window will depend on the observable in 
question. GW fermions help a lot since only in this case we have a variant 
of chiral symmetry at finite lattice spacings while strictly speaking in other 
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Figure 2. Chiral extrapolation of the nucleon mass, after phenomenological corrections 
for finite lattice spacing and volume effects. The data are from the CP-PACS'^^ (solid 
circles) and UKQCD Collaborations^* (open circles). 



formulations (without doublers) a sensitive comparison can only be per- 
formed after extrapolation of lattice results to the continuum limit. While 
even 400 MeV <^ ^xSB (modulo the ambiguity of vs. "<") such a pion 
is still doomed to "see" some of the internal structure of the proton, with 
an inverse charge radius of about 250 MeV: it is doubtful that the quark 
and gluon nature of QCD can completely be ignored with such a "hard" 
pion probe. With broadening baryonic wave functions such issues are likely 
to be more critical in QCD with light sea quarks than they already are in 
the quenched approximation. 

Naively, hadron masses are a polynomial in the quark mass, mg oc m^. 
However, pion loops give rise to non analytic (NA) functional dependence 
on the quark mass. For instance the nucleon mass is given by, 

mAr(m^) = TOjv(O) + C2m^ + csTO^ H . (1) 

The coefficients of the NA terms are related to phenomenological low en- 
ergy constants. For instance C3 = — 3.g^/(327rF^). In the quenched approx- 
imation, which neglects disconnected fermion lines, the leading NA (LNA) 
term is not proportional to anymore but to InTOjr, known as a chiral 
log. Only very recently convincing signs of such logarithmic behaviour have 
been found, albeit not yet with the expected coefficients. 

Recently, corrections to standard LNA terms due to massive pions (re- 
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suiting in a different pion cloud) have been modelled, in particular by the 
Adelaide group. Unlike the chiral Lagrangian itself, such models cannot 
be derived in any limit directly from QCD, and the predictive power has to 
be digested with some caution. Figure 2 reveals that wc have not yet made 
contact with chiral perturbation theory (dotted curve, C3 fixed). Nonethe- 
less the data are well fitted by Eq. (1), with C3 as a free fit parameter 
(dashed curve), and the extrapolated value has dangerously tiny statisti- 
cal errors. The Adelaide cloudy bag model (CBM) suggests a different 
functional form that screens the LNA term at pion masses close to a bag 
parameter, A. With a fixed phenomenological value of A the solid curve 
is obtained while a three parameter fit yields the dashed-dotted curve. By 
treating the difference as a systematic uncertainty such models allow us to 
arrive at more realistic error estimates. Alternatives to the CBM include 
NJL based^^ as well as Bethe-Salpeter based^^ ideas. 

The quenched approximation is sick in several ways. I already men- 
tioned that many unstable particles become stable. It is not always entirely 
clear when this "bug" is a virtue and in which cases not. A more serious 
flaw is that unitarity is violated since quarks are not included in a consis- 
tent way. This defect, which becomes visible in particular at small quark 
masses and in the scalar sector, also reflects onto and can be understood 
in terms of the chiral expansion. For instance chiral logs originate from 
this source. Another related problem is the absence of the axial anomaly 
since quarks do not feed back onto the glue. This leads to the spontaneous 
breaking of a greater flavour symmetry group to start with and results in 

instead of njp — 1 pions: the rj' becomes heavy only at small sea quark 
masses, with possible implications onto the phenomenology of excited 
baryon states: the mass splitting between the nucleon {N = N^^^ ) and 
the N* = N^/'^ for instance is a consequence of chiral symmetry breaking 
and its size is related to transitions N* N + tt ^ TV*. The number 
of pions should naturally be expected to seriously affect such observables. 
(Note that also in the quenched approximation there is a pion cloud since 
valence quarks can travel forwards and backwards and forwards again in 
time. Pion exchange is possible too, as valence quarks are swapped.) 

To this end it has recently been suggested to model quenched LNA 
terms and to subtract these from the corresponding lattice simulation re- 
sults, replacing them by the cxpectcxi 'u,?7,-guenc/ied counterparts.-'^* Needless 
to say that many assumptions arc required, in particular relations between 
low energy parameters appearing in quenched and un-quenched chiral La- 
grangians and in the CBM. However, in the absence of lattice results with 
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light sea quarks, playing around with such models provides us with valuable 
rough ideas of possible sizes of quenching effects. On the other hand, once 
quality un-quenched data become available this will positively feed back 
onto such models, providing insights impossible to gain from experiment 
alone where quark masses and the number of active flavours are fixed. 

Determinations of nucleon parton distributions seem to suffer from ex- 
tremely large uncertainties in the chiral extrapolation; in particular this 
is so for the leading moments. In a recent study of the LHPC/SESAM 
Collaborations^^ of quark distributions including sea quarks almost per- 
fect agreement was found with quenched results, obtained both by them 
and by QCDSF,^" down to « 600 McV. However, a naive polynomial 
extrapolation of {x)u-d was found to overestimate the phenomenological 
value by as much as a factor 5/3. Such a factor can indeed be accounted 
for by modelling pion mass effects on the LNA contribution in a sensible 
way,^^ indicating that the quark mass is the dominant (and a huge) source 
of uncertainty. 

4. Spectrum 

In the past two years there has been a vibrant flow of lattice publications 
on the spectrum of excited nucleon states. All but one^^ have been based 
on the quenched approximation so far. Only the LHPC/UKQCD/QCDSF 
Collaborations,^'* using the Wilson-Clover action, attempted a continuum 
limit extrapolation. Other strategies were implementations of improved ac- 
tions like the D234 action by Lee et al?^ or the FLIC action by Melnitchouk 
et alP'^ . The Rikcn-BNL group employed domain wall fermions to improve 
the chiral properties. 

One feature that all these articles share is a level inversion between the 
positive parity Roper A^'(1440) resonance and the negative parity A''* (1535), 
relative to experiment. It was not clear whether this was a quenching effect, 
hence shared with many models, most of which are essentially quenched 
too. But maybe the resonance observed in nature did something funny, 
to which the quark-only operator used on the lattice to create the Roper 
did not couple well. A pessimist might also argue that such resonances are 
quite unstable which could very well affect their position in ways that are 
not easily foreseeable, prior to doing the real thing. 

However, new results obtained with GW fermions by the Kcntucky- 
Washingon group^'' at lighter pion masses, > 180 MeV, than ever re- 
alised before indicate an exciting behaviour: as the pion becomes lighter 
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than about 550 MeV the slope of the Roper data bends downwards, head- 
ing straight in the direction of the experimental value (cf. Figure 3). While 
at the lightest quark mass LrriT^ Ki 2.7 and finite size effects are likely, ev- 
erything with rn^ > 0.1 GeV^ is pretty safe. There is overlap with other 
lattice studies^^'^^'^^ in the region m| > 0.25 GeV^ and one might wonder 
why the onset of this behaviour has not been detected before. Inspection 
shows, however, that within the statistical errors an effect of the expected 
size would not have appeared significant in these cases. 

Certainly, once the independence of the Kentucky results from Bayesian 
fitting priors and lattice volume has been firmly established, these are very 
exciting news that tell model builders how transitions between heavy and 
light quark regions can look like: a prime example for the cross fertilisa- 
tion between lattice QCD, models and experiment in this thriving area of 
research. The Kentucky group has investigated many additional baryonic 
channels and it is astonishing how good quenched QCD seems to explain 
the positions of these strongly decaying resonances, once the quark masses 
become sufficiently light, another unsolved miracle of nature. 
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Figure 3. Chiral extrapolation of the A', N' and A^* masses. 
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5. Structure 

Exciting results on the transition matrix element for electromagnetic N to 
A transitions have been reported. The relevant operator is usually decom- 
posed into three Sachs type form factors, accompanying magnetic dipole 
transitions, Gmi{<i'^) as well as electric and Coulomb/scalar quadrupole 
transitions, Ge2{q^) and Gc2((?^)- The quadrupole transitions ratios 
Rem = Gei/Gmi and Rsm = GciIGmx are related to the question of 
nucleon deformation.^^ Recent experimental results, in particular from the 
CLAS^^ and OOPS Collaborations^o indicate the ratios Rem ~ -0.02 
and Rsm ~ —0.07 around « 0.53, which is the lowest momentum 
transfer that could be realised in a recent lattice study with two flavours 
of sea quarks, not much lighter than the strange. "^'^ These first quantita- 
tive lattice results seem to get the sign right and point towards values 
Rem = —0.03(1) and Rsm = —0.03(2), consistent with experiment. In- 
terestingly, the quenched reference result Rem = —0.009(8) is compatible 
with zero, indicating that a more realistic pion cloud might be essential 
for the effect. Obviously, lattice simulations are still miles away from com- 
puting, say, octupole form factors for transitions between orbitally excited 
nucleons. 

In another lattice investigation^^ the electric charge and matter distri- 
butions of TT, p, N and A have been studied at rest. One can easily calculate 
the charge radius from such a distribution and the Fourier transform is re- 
lated to form factors at zero momentum transfer. While the (cut-off and 
scheme dependent) matter distributions of all investigated hadrons were 
found to be remarkably similar there are differences in the charge distribu- 
tions. In the quenched approximation the pion is only approximately half 
as wide as the A, with A'' and p almost of the same size as the A. Even 
with a pion mass of almost 600 MeV un-quenching significantly broadened 
the latter three distributions, and in particular the A, while the tt remains 
completely unaffected. The quenched proton charge radius comes out to 
be Tp = 0.59(4) fm, after polynomial extrapolation to the physical limit. 
This has to be compared with the value Vp » 0.81 fm from dipole fits to 
experimentally determined electromagnetic form factors. The above direct 
lattice result is in agreement with rp « 0.6 fm obtained indirectly from a 
dipole fit to electric and magnetic proton form factors calculated on the 
lattice by the QCDSF Collaboration.^^ However, from recent JLAB Hall A 
results, "^^ showing an energy dependence of the ratio Ge/G„i, we know that 
the dipole ansatz is not the whole truth, which stresses the importance of 
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determining quantities like charge radii in more than one way. The same 
QCDSF work suggests rp « 0.7 fm for two flavours of sea quarks, closer to 
the phenomenological value. 

For the first time deformations of the hadron wave functions have been 
investigated.'^^ Obviously, a spectroscopic quadrupole moment becomes 
only possible for an angular momentum of at least one, such that devi- 
ations from spherical symmetry are ruled out for the nucleon. However, 
the p was reported to be somewhat prolate with the A being slightly oblate 
and the deformations appeared to increase when including sea quarks. In 
neither of the cases was the deformation of the A statistically significant 
though. 

It is not clear whether QCD forces are stringy with a three body po- 
tential that depends on the length of the shortest path connecting the 
three quarks to one point (Mercedes or Y law) or if they merely can be 
understood in terms of sums of two body interactions (A law)."''"' Leaving 
aside the question whether an instantaneous Born-Oppenheimer interaction 
picture is appropriate for a situation with light quarks, the nucleon wave 
fimction seems to be better described by a Schrodinger equation with a 
potential that depends on the "A-distance" than by the Y configuration.^^ 
This finding is also supported by recent lattice calculations of SU{^) and 
SU{A) gauge theory as well as of the Potts model. The results indi- 
cate that for interquark distances of physical relevance, i.e. r < 1 fm, the 
A-picture yields a very accurate description while, long after the "string" 
of QCD with sea quarks is broken, the Y picture is eventually approached. 

Recently, spin and transversity of the quark contributions to the A hy- 
peron have been calculated by QCDSF. '^^ The fraction A^a of the spin 
carried by quarks of fiavour q is given by the forward matrix element of the 
axial vector current, (A(p', s)|g7^75g|A(p, s)) = 2 (p' — p)^Aq\, i.e. by the 
lowest moment of the hclicity distribution. Assuming SU (3) flavour sym- 
metry one can express the AgAS in terms of proton spin fractions Ag^. Con- 
verting the lattice results into the MS scheme at /x = 2 GeV yields an axial 
charge = Aup — Adp that falls short of the experimental value ~ 1-26 
by 10-20 %, due to quenching and/or due to uncertainties in the chiral ex- 
trapolation, such that we cannot expect the result Aua = Ac^a = —0.02(4) 
to be consistent with the real world either. 

The lattice data on proton and A spin fractions turn out to be in ex- 
cellent agreement with SU (3) flavour symmetry, an observation that might 
be independent of sea quarks and certainly consistent with the success 
in understanding semileptonic decay rates of the A. Combining experi- 
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mental data obtained for the proton helicity with the systematic error on 
SU{3)f symmetry suggested by this lattice study, yields the predictions 
AuA = AdA = -0.17(3) and Asa = 0.63(3), again in the MS scheme at 
H = 2 GeV. 

6. Outlook 

Experimental efforts are paralleled by lattice studies of baryon spectrum 
and structure. The main limitation that has to be overcome in the near 
future is the use of unrcalistically heavy quarks in quenched simulations, 
as well as in simulations with sea quarks. Recent theoretical developments 
and advances in computer technology mean that optimism is justified. 
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